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ABSTRACT 

This work deals with approximated performance 
evaluation for airborne navigation systems composed of 
GPS and low cost inertial sensors. Usually, 
performance evaluation requires to design a 
multivariable Kalman filter and to perform intensive 
numerical studies. This work follows [1] and proposes 
to consider an INS/GPS integration algorithm as a 
fusion of two processes. The first one is called filtering 
and deals with the sensor error model. It is aimed to 
calculate the high order position derivatives (velocity, 
acceleration, jerk). The second process is named 
separation and its goal is to calculate, from filtering 
results, all navigation and sensor errors. The most 
important feature of the above partition is the ability to 
analyze independently sensor error models and flight 
trajectories. 

This paper analyzes the accuracy achieved by 
the calculations based on filtering-separation partition. 
First the underlying assumptions are precisely defined. 
The errors introduced by those assumptions, for 
different trajectories and sensor models, are analyzed. 
Algorithms for filtering and separation are discussed. 
Finally, leveling errors calculated by the combined 
filtering-separation algorithms are shown to be very 
close to the ones calculated by the standard, complete 
Kalman filter. 

1 INTRODUCTION 

An integrated Global Positioning System (GPS) 
receiver with an Inertial Navigation System (INS) is 
widely used as a precise, high bandwidth, robust, 
position, velocity, and attitude data source. See 16], [7] 
for a general reA^iew and discussion about benefits of 
INS-GPS systems. A standard technique to integrate 
INS with GPS is to design an appropriate Kalman 
filter. The integrated system is highly coupled, it is 
time, trajectory, and sensor error model dependent. 
Therefore, mainly because of the need to carry out 
intensive numerical studies, the task of designing a 
Kalman filter is time consuming. Consequently, a tool 
to evaluate approximately, but easily and fast, the 
system performance is highly desirable. This paper 
makes an attempt to create a general framework, in 
which the influence of the sensor error model and the 
flight trajectory are analyzed easily and independently. 

The idea is to divide the entire estimation 
algorithm into two processes: filtering and separation. 



Filtering is aimed to estimate quantities like position, 
velocity, acceleration, and jerk. They are described by a 
time-invariant system. This system takes into account 
the statistical sensor error model. The second process, 
assumes as known the results of filtering, mainly 
acceleration and jerk, and calculates all unknown 
navigation and sensor errors, for example attitude, gyro 
drifts, or accelerometer biases. This process is called 
separation and it deals with the flight trajectory. 

The goal of this work is to quantify how well, 
for different trajectories and sensor error models, this 
partitioning works. The paper is organized as follows. 
First the standard strapdown navigation error model is 
presented. Some model simplifications are proposed 
and justified for medium (tactical) or low (automotive) 
quality of inertial sensors. The next proposition is to 
consider the flight trajectory as consists of level straight 
legs, with short turns connecting consecutive segments. 
Our performance criterion is attitude estimation 
accuracy. We conclude that for a general trajectory, the 
partition into equivalent straight segments is related to 
properly identified time constants. Those time 
constants are associated with acceleration and can be 
calculated by filtering problem analysis. Finally the 
filtering and separation algorithm are discusses, their 
combined results are successfully compared with the 
standard Kalman filter calculations. 



2. MATHEMATICAL BACKGROUND 

The standard error model for a strapdown system is as 



follows: 

r~w r xr + v (1) 

v = w v xv + fxyr + 8g + C BL B (2) 

\y = w ¥ XY-C BL D (3) 



where r^lff are the position, velocity, and attitude 

error vectors, respectively. f y Sg, C LB are the specific 

force, gravity error, and transition matrix from sensor 
(body) to navigation (LLLN -local north local level) 
frame. B,D arc the bias and drift error vectors, 
respectively. 

The angular velocities are defined by: 
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w ,, 



(6) 



w_ LE y w m are the rotation vector from navigation 

frame to earth fixed frame and the rotation vector from 
the earth fixed frame to the inertia! frame, respectively. 
See [5] for more details. 

For medium (tactical) and low (automotive) 
quality inertia! sensors the time of concern is relatively 
short, bias, drift, and attitude errors are relatively high, 
and, as it turns out one can neglect the terms 
containing rotation vectors and the gravity error. Doing 
so, we obtain the following approximated navigation 
error model. 

t = V (7) 

v = fxW + C LB B (8) 



(9) 



Let us discuss the dynamic model of drift and 
bias. For proper GPS/INS integration, it is very 
important to model correctly variations in bias and 
drift. The standard approach is to define two 
components, random constant for constant drift or bias, 
and Markov process for drift (bias) variations. Let 

D k \B k denote constant drift and bias at discrete 

time t k9 similarly D k m , B k m denote discrete time 

Markov drift and bias. Then the discrete time system 
equations are written as follows: 

>U=r k +v k dt (io) 
•■v k +F k yr k +C k B k '+C k B k m +w k y (li) 



D M "=D k 



D kM -a D k +w k 



B t J =b: BJ v =a B B™ + w b B 



(12) 

where F k is the integral, from t k till t k+l , of the skew 
symmetric matrix corresponding to the specific force 
/. Similarly, C T A . is the integral of the transition 

matrix C LB . Noise processes w k \w k v are due to 
accelerometer and gyro random walk, respectively. The 
time increment is denoted by dt . 

In order to define precisely the filtering process, let us 
define the following new variables. 



A k :=^(F k¥k +c t B;) 

at 
at" 



(14) 

(15) 
(16) 



■h m =~F kM C k D k m (17) 
at 

For constant flight segments (f r k+2 = F M = F k ) 

(Pk+i = ) liew variables behave like constant 
acceleration, Markov acceleration, constant jerk and 
Markov jerk. Namely, it turns out that they satisfy the 
following equations: 



r M =r k + v k dt 

= v k +A k G dt + A k m dt + w k v 

A k+ ;=A;+j;dt+j k m dt+w k J 

m s^B a rn 



A m B A m . > 



/ c = J c 



M =a u J k +w k 



(18) 
(19) 
(20) 
(21) 
(22) 
(23) 



where w k , w k ,w k A ,w h , m are corresponding 
components of process noise. 

Observe that Equations (18-23) are trajectory 
invariant, moreover they describe time constant system. 
On the other hand they cover all sensor error model 
parameters. Equations (14-23) define the division of 
estimation algorithm into two processes. The filtering 
process is defined as an integration of Equations 
(18-23) with the GPS error model. The separation 
process is defined as a solution of Equations (14-17) 
with left hand sides assumed, as filtering process 
results, to be known. Note that the flight trajectory 

data, described by F k and C k appear only in 

separation problem set-up. All data about the sensor 
error model, process noise and Markov process 
parameters are included only in filtering problem 
description. 

The above approach is abased on two 
assumptions: 

• The navigation error model is described by 
Equations (7-9) 

• The flight trajectory can be described as consists of 
constant segments connected by short-time turns. 

In the fol lowing section, we describe numerical studies 
to justify the above assumptions. 

3. PRELIMINARY NUMERICAL ANALYSIS 

To carry out numerical studies we need to 
define precisely the error model and flight trajectories. 
We consider a parametric family of IMU's (Inertial 
Measurement Units) named Ql, Q10, Q100. The name 
convention is that Ql is a prototype of a standard 
tactical quality IMU (for example, HonyweH's FIG 1700, 
Litton's LN200, or Kearfotfs T16B), Q100 is a 
prototype of automotive quality sensors. Q10 is a 
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theoretical model of a unit roughly 10 times worst than 
Ql. Table 1 describes their detailed error models. All 
Markov processes defined in Table 1 have 300 seconds 
correlation time. 



Table 1 IMlTs error model 



Error 


Ql 


Q10 


Q100 


gyro constant drift 


1 deg/h 


10 deg/h 


100 deg/h 


gyro Markov drift 


0.2 deg/h 


2 deg/h 


20 deg/h 


gyro random walk 
deg/sqrt(h) 


0.1 


0.3 


1.0 


accelerometer 
constant bias 


lmg 


3mg 


10 mg ! 


accelerometer 
Markov bias 


0.2 mg 


0.6 
mg 


2 

mg 


accelerometer 
random walk 
m/sec/sqrt(h) 


0.03 


0.1 


0.3 



GPS error model is given in Table 2. From 
many possible configurations, without effecting the 
final conclusions, SA off and position measurements 

are assumed 

Table 2 GPS error model 



Error 


Variance 


Correlation 
time 


Horizontal position 
white 


2.8 m 




Horizontal position 
Markov 


10 m 


20 sec. 


Vertical position 
white - 


3.5 m. 




Vertical position 
Markov 


■13 m. 


20 sec 



The flight trajectories, analyzed in this paper, are level 
and consist of straight segments and turns. They have 
three parameters. The angle a is the change of the 

heading during turns, t m \s the time between two 
consecutive turns, and t d is the turn duration. 

Having defined all error models (EvlU and 
GPS) it is straightforward to define the appropriate 
Kalman filter. This filter calculates the integrated 
system optimal performance. We will focus on the 
standard deviation of attitude errors, leveling yr x and 
heading \j/ s . Tables 3.1 and 3.2 present results for half 
an hour simulation with time between turns 
/ =300 seconds. 



Table 3.1 Leveling error for a = 15 



w 

1 X 


Ql 


Q10 


Q100 


mrad 














u 


Prec. 


Appr 


Prec. 


Appr 


Prec. 


Appr 


1 


0.58 


0.70 


1.71 


2.08 


5.59 


6.83 


5 


0.58 


0.70 


1.71 


2.09 


5.60 


6.85 


10 


0.58 


0.70 


1.71 


2.09 


5.66 


6.91 


20 


0.58 


0.71 


1.73 


2.11 


5.85 


7.13 


50 


0.59 


0.72 


1.83 


2.23 


6.96 


8.23 


100 


0.64 


0.77 


2.07 


2.47 


8.64 


9.49 



Table 3.2 Heading error y/ z for a = 15 



mrad 


Ql 


Q10 


Q100 


u 


Prec 


Appr 


Prec 


Appr 


Prec 


Appr 


1 


7.8 


8.0 


23.1 


23.4 


51 


51 


5 


7.8 


8.1 


23.5 


23.8 


53 


53 


10 


7.9 


8.1 


24.1 


24.4 


55 


56 


20 


8.1 


8.3 


25.6 


26.1 


62 


62 


50 


8.9 


9.2 


33.0 


33.8 


102 


103 


100 


10.6 


11.2 


49.6 


51.7 


237 


242 



The column "Prec." corresponds to the 
simulation based on the precise error model (equations 
(1), (2), (3)) while the column "Appr." is related to 
approximated error model (equations (7), (8) (9)). 
From the comparison between corresponding columns, 
one can see that using the approximated strapdown 
error model introduces differences less than 20 %. This 
observation has twofold consequences: 

• It justifies the approximated model 

• It sets the accuracy level that should be expected in 
the subsequent analysis 

It is interesting to consider the influence of 
turn duration t d on the performance. It is clear from 
Tables 3.1 and 3.2 that for low values, 1-10 seconds, 
t d does not influence the results. In other words, two 

different trajectories, with the same heading changes 
and the same time between turns, but with different 
turn duration time, yield the same performance. This is 
key observation for the idea to divide the estimation 
algorithms into two parts: time invariant filtering and 
trajectory dependent separation. Since for separation a 
flight trajectory with constant segments is required, the 
idea is to replace every turn with an equivalent one 
with very short duration time. The next step is to 
neglect measurements during those very short turns, 
namely to consider estimation, as in separation set-up, 
only during constant segments 

In spite of the fact that for low t d the above 
idea works very well, observe that for longer duration 
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times (> 20 seconds) there are some changes in 
performance. Therefore, it is very important to find 
maximal duration time that is still equivalent to a very 
short duration time, say f d ~\ second. Figure 1 
describes the normalized leveling error with respect to 
turn duration time t d . The heading change is 45 

degree. It is evident that even the normalized behavior 
depends on sensor quality 

Table 4 compares two different numbers. The 
first one is the maximal duration time to preserve 
leveling accuracy within the limits of 10 percents. To 
explain the second number, consider the time-invariant 
model, described in equations (18-22) integrated with 
GPS, as defined in Table 2. The resulting Kalman filter 
achieves steady state, and then the integrated system is 
time-invariant as well. The second row in Table 4 
describes the time constants for acceleration in this 
steady state system. 

normalized psix 

1-7 j 1 1 1 1 — i 1 1 1 1 1 



1.6 -• 




0 10 Z0 30 40 50 60 70 80 90 100 



Figure 1 Normalized \jf x 



Table 4 Characteristic times comparison 





Ql 


Q10 


Q100 


t d lor 10% accuracy 


> 100 


80 


40 


Ti me constant for 
acceleration in 
time-invariant svstem 


121 


75 


36 



The conclusion is apparent: as long as we deal 
with duration times less than time constants of the 
filtering process, we can reduce the duration time 
without changing significantly the performance. This 
conclusion has been further verified for a wide range of 
problem parameters. Now we can state the algorithm 
for approximate performance evaluation. 



• Analyze the performance and time constants of 
filtering problem with respect to sensor error 
model. 

• Divide the actual trajectory into segments no 
longer that the acceleration time constant as 
calculated for filtering problem. Preserve the 
heading change for every segment. 

• Combine the results of filtering problem for 
different segments (separation process) 

4 FILTERING AND SEPARATION 

The solution of filtering problem is straightforward, it 
is carried out by Kalman filter which integrates 
dynamics described by Equations (18-23) with GPS 
error model specified in Table 2. The remaining 
separation problem is to solve Equations (14-17) 
assuming that: 

• the flight trajectory consists of L constant 
segments, 

• estimates and error covariance of accelerations 
and jerks [A k ° , A k m , J k c , J k m ) are given. 

Since flight trajectory segments are separated by time 
greater than filtering time constant, we can assume that 
acceleration and jerk estimates for different segments 
are statistically independent. This key observation 
provides the means to define a simple Kalman filter 
with attitude, biases, and drifts as states and 
acceleration and jerk as measurements. This approach 
yields accurate results, but further problem 
simplification is possible. 

Consider again Equation (14), bias B k ° is 

constant but the attitude \ff k , due to gyro drifts and 

gyro random walk, is time dependent. On the other 

hand, when acceleration estimate A k was calculated, 

drift influences have been taken into account. 
Therefore, perhaps, one can consider a mean value 
version of leveling error. Namely, let us assume: 

4 e +v/=^(F^+c t zr) k=u 2 ,...,i L 

where A k 7 F k ,C k ,dt are given. The additional 

measurement noise v k 8 is due to gyro random walk. 

The covariance matrix of the measurement errors R 
consists of a I , due to acceleration errors, and W due 
to gyro random walk. The problem has two unknown 
constant vectors \jf,B c . This is a linear constant 

parameter estimation and least square approach is the 
proper solution. Lot us define 
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0 10 — fstaricleircJ cJeviatlon of psi x [mrad] 
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KF F-S 


V\ 




PS — precise 




40 so 

alpha Ideg] 



Figure 2 Leveling error with respect to heading change CC - as calculated by different methods 



T = 



(24) 



It is well known that the covariance matrix of the 
unknown parameters is given by: 



p=fc m T R-*T M y 



(25). 

Figure 2 describes the standard deviation of 
leveling error calculated by three means: 

• LS F-S (least square filtering-separation) as 
described by Equations (24-25) 

• KF F-S (Kalman filter filtering-separation)- 
Kalman filter solution of separation problem. 

• PS (precise) - calculated by the complete 
simulation. 

The first two results have been modified with respect to 
the following observations. The steady state leveling 
error can not be greater then the initial accelerometer 
bias error (multiplied by the gravity). The leveling 
error can not be less then the acceleration error 
(multiplied by tire gravity) as calculated by the filtering 
problem 

From Figure 2 one can conclude that the 
attitude accuracy, calculated by the filtering- separation 
scheme with a low dimensional Kalman filter, is very 
close to the one obtained by the entire system 
simulation. The above is true for all tested trajectories, 
those with large a , namely good estimability (see [3], 
[4] for an interesting discussion about observability of 
piece-wise constant systems) and those with small 
heading changes (poor estimability). Figure 2 depicts 
the results for Q10, the plots for Ql and Q100 are very 



similar. The results of least square are very good for 
large heading changes, but are slightly worst for low 
estimability. Indeed, the least square is not able to take 
into account initial covariance, which, in this case, is 
an important factor. 

5 CONCLUSION 

It has been shown that the idea to divide the estimation 
process into filtering and separation provides simple 
and quite precise meanings for system performance 
evaluation. The most important feature of this approach 
is the independent treatment of the sensor error model 
and the flight trajectory. 
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